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We study the dynamics of a non-minimally coupled scalar field cosmology with a potential function. We 
use the framework of dynamical systems theory to investigate all evolutional paths admissible for all initial 
conditions. Additionally, we assume the presence of barotropic matter and show that the dynamics can be 
formulated in terms of an autonomous dynamical system. We have found fixed points corresponding to 
three main stages of the evolution of the universe, namely, radiation, matter and quintessence domination 
epochs. Using the linearization of the dynamical systems in the vicinity of the critical points we explicitly 
obtain formulas determining the effective equation of state parameter for the universe at different epochs. 
In our approach the form of w(z) parametrisation is derived directly from the dynamical equations rather 
than postulated a priori. 
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There are principally two approaches in searching for the form of the equation of state (EoS) for the 
current Universe in the accelerating phase of expansion. In the first approach the form of the coefficient 
Wx(z) of the equation of state for dark energy is assumed at the very beginning, usually in the simple 
linear form with respect to the redshift z or the scale factor a 0][2] ■ In the second approach the form of the 
EoS parameterisation is derived from the exact dynamics of the underlying theoretical model SJ7). In the 
present paper we realise this idea for the Friedmann-Robertson- Walker (FRW) models filled with a scalar 
field non-minimally coupled to gravity. 

In our study of evolutional scenarios it is natural to use the framework of dynamical system theory 
because it offers the possibility of investigations of all solutions admissible for all initial conditions. We 
investigate fixed points of a dynamical system and then linearise the system around them to find exact forms 
of the EoS parameterisations. Therefore we derive w{z) parameterisations directly from the underlying 
dynamics of the model. 

We assume the spatially flat FRW universe filled with a non-minimally coupled scalar field and barotropic 
fluid with a general equation of state parameter w m , and the action is 



\r - e (g^d^d^ + tiRtf) - 2U{4>) 



S m , (1) 



where k 2 = 8ttG, e = +1,-1 corresponds to canonical and phantom scalar field, respectively, and the 
metric signature is (— , +. +, +). S m is the action for the barotropic matter part. 

We can obtain the dynamical equation for the scalar field from the variation SS/S<p = 0, and energy 
constraint from the variation 5S/ Sg = : 

4> + 3Hcj) + £,R(t ) + eU'(<t)) = 0, £ = e\ft + e3£flV + e^H{(j) 2 )+ U{<t>) +p m - ^H 2 . (2) 

2 K z 



* E-mail: hrycyna@kul.lublin.pl 
** E-mail: uoszydlo@cyf-kr.edu.pl 



Copyright line will be provided by the publisher 



2 



O. Hrycyna and M. Szydlowski: Three steps to accelerated expansion 



Table 1 The location and eigenvalues of the critical points in twister quintessence scenario 



WcS 


location 


eigenvalues 


1 
3 

w m 

-1 


x* 2 = 0,y* 2 =0,\* 2 = 

x* 3 = o,(y* 3 ) 2 = i,K = o 


h = 

h,3 = "1(1-1 

h,a = -i 


= -66 z 2 = i2f, i 3 = * 

^3 ± ^9 + e2a - 48^ 


5£(1 - 3ti) m ) 
E^).fa = |(l + t« m ) 
,h = -3{l + w m ) 



Then the conservation conditions read 

3 • n 2 

—H 2 = P4. + p m , H = — — [(pj + p^) + p m (l + w m )] , (3) 
where the energy density and the pressure of the scalar field are 

P^ = s\i> 2 + U(<j>) + eZtH 2 4> 2 + s^H{4> 2 ), (4) 

P<t> = 4 (1 ~ 4e) ^ ~ U W + ^H(d> 2 ) -e2^(l - 6£)H</> 2 - 63^(1 - 8^)H 2 <f> 2 + 2^U'(<f>) ■ (5) 
In what follows we introduce the energy phase space variables 



k4> K\/U{6) 



V3H 



(6) 



which are suggested by the conservation condition jfj^Pj, + wjjzp m = ^0 + = 1. 
The acceleration equation can be written in the form 

k 2 3 

H = — (peff + Peff) = -^H 2 (l + fflcff) , 

where the effective equation of state parameter reads (a = — uj^) & d^f ^ ) 
1 



(7) 



l"cff 



1 - e6f (1 - 6£)z 2 



U{4>) d<t> J 

1 + e(l - 60(1 - w m )x 2 + e2f(l - 3w m )(x + z) 2 + 



+ (1 + w m )(l - y 2 ) - £2^(1 - Q£)z 2 - 2^\y 2 z 
The dynamical system of the model under consideration takes the form JUS! 

x' = -3.t - 12£z + e-Xy 2 1 - e6£,z(x + z) + e6£(l - &£)xz 2 + 
+ l(x + 6£z) 



(8) 



3 

+ 2 lJ 



z = x 



£ (1 - 6^(1 - w m ).T 2 + £2£(1 - 3w m )(x + z f + (1 + w m ){\ - y 2 ) 
1 - e6£(l - 6C)z 2 + 6&A j - el2f (1 - 6£)yz 2 + 
e(l - 60(1 - w m )x 2 + £2^(1 - 3w m )(a; + z) 2 + (1 + w m )(l - y 2 ) 
1 - e60l - 60-z 2 



(9) 

(10) 
(11) 



A' = -A 2 (T -l)x 1 - £601 - 6£)z 



(12) 

Copyright line will be provided by the publisher 



adp header will be provided by the publisher 



3 



where T = and a prime denotes differentiation with respect to time r defined as — 



1 



e6£(l — 6£)z 2 dl „ n - In the rest of the paper we will assume the function T in the form T(X) = 1 — -p-. 
In our approach the exact formulas for w{z) parameterisations can be derived directly from the linearised 
solutions of the dynamical system in the vicinity of the critical points representing different epochs in the 
evolution of the universe (Table[T]i. In the special case of conformally coupled (£ = 1/6) canonical (e = 1) 
scalar field and dust matter (w m = 0), they are as follows: 

1) At the radiation domination epoch, 



where A x = xf - ± (x^ - A*), B x = xf + ± - A*), (A?) 2 = a 2 and z R is the redshift of the 
radiation domination epoch. 

2) At the matter domination epoch, 

(i) A^ (z) A^ 

where ^2 = ^2 + ~2~> -»2 = #2 H — — and zm is the redshift of the matter domination epoch. 

3) In the de Sitter state, we have two parameterisations whose type depends on the value of the parameter 
a characterising the shape of the function V. 

— For a > —1/2, we have a linear approach to the de Sitter state 

«&(*)!*=* = -1 + g£ - a + ^) ^3(1 + ^) 3+v3 + 

+ (|A 3 B 3 - 2yl {yf 2/3*)) (1 + *) 3 + A (l - a - Va) B 3 2 (l + zf~^- 

^Ajy* (yW _ y *) (1 + z) a + 7A + (g^tf - ^ + tf) (y? - tf) (1 + zf- 

- H 1 ^ - vi) (1 + -) 6 ^- 

- ^ - y 3 *) 2 (4,(1 + + B 3 (i + z)-^) 2 (1 + z) 9 (15) 
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Fig. 1 The evolution of w c g given by the relation ^ for the non- 
minimally coupled canonical scalar field e — +1 and the positive 
coupling constant f. The existence of a short time interval during 
which w e g — | is the effect of the nonzero coupling constant £. 

where A 3 = xf + ^(3 - \/A)\$\ B 3 = x£> + ^(3 + and A = 1 + 2a > 0. 

— For a < —1/2, we have a damped oscillatory approach to the de Sitter state 



«& = -i - 2»5 [vf - vi) (i + ^) 3 + 



where C; 




These parameterisations depend on three values Xn \ Vn \ An of the initial conditions for the linearised 
solutions at different epochs and the parameter a describing the shape of the function T(A). The common 
parameter a is present in all parameterisations and can be estimated from the observational data. 

In this short note, we presented the possibility of extracting the equation of state parameterisations of 
dynamical dark energy directly from the dynamics of the underlying theoretical model. It is interesting that 
the obtained formulas for w(z) are all that is needed for the realistic cosmological model, i.e., the radiation 
epoch required by the nucleosynthesis, matter domination phase and final acceleration epoch (see Fig.[T|). 
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